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ABSTRACT 



J^j ■ We present high resolution (R ~ 44,000) UVES spectra of the eclipsing pre- 

main sequence star KH 15D covering the wavelength range 4780 to 6810 A ob- 
tained at three phases: out of eclipse, near minimum light and during egress. The 
system evidently acts like a natural coronagraph, enhancing the contrast relative 
to the continuum of hydrogen and forbidden emission lines during eclipse. At 
maximum light the Ha equivalent width was ~2 A and the profile showed broad 
wings and a deep central absorption. During egress the equivalent width was 
much higher (~70 A) and the broad wings, which extend to ± 300 km/s, were 
prominent. During eclipse totality the equivalent width was less than during 
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egress (~40 A) and the high velocity wings were much weaker. B./3 showed a 
somewhat different behavior, revealing only the blue-shifted portion of the high 
velocity component during eclipse and egress. [01] AA6300, 6363 lines are easily 
seen both out of eclipse and when the photosphere is obscured and exhibit lit- 
tle or no flux variation with eclipse phase. Our interpretation is that KH 15D, 
although clearly a weak-line T Tauri star by the usual criteria, is still accreting 
matter from a circumstellar disk, and has a well-collimated bipolar jet. As the 
knife-edge of the occulting matter passes across the close stellar environment it is 
evidently revealing structure in the magnetosphere of this pre-main sequence star 
with unprecedented spatial resolution. We also show that there is only a small, 
perhaps marginally significant, change in the velocity of the K7 star between the 
maximum light and egress phases probed here. 

Subject headings: starsdndividual (KH 15D) — accretion — line:profiles 

1. Introduction 

KH 15D is a unique system in which a pre-main sequence star is periodically occulted 
by extended, non-luminous matter, presumably part of a circumstellar disk (Herbst et al. 
2002). Every 48.36 days the star fades over 2-3 days by ~3.5 mag and remains near minimum 
light for ~20 days. The eclipse duration has been increasing with time, by ~l-2 days/yr. 
Hamilton et al. (2001) obtained low resolution spectra in and out of eclipse and concluded 
that the star was a K7 weak-line T Tauri star (WTTS). The WTTS classification, based 
on the equivalent width (EW) of the Ha line out of eclipse, is supported by the lack of any 
IR 2 or UV excess emission and the absence of substantial photometric variability outside of 
eclipse. Based on its membership in NGC 2264 (Sung et al. 1997), the estimated mass of 
the K7 star is 0.5-1 M and its age is 2-4 Myr. 

In an attempt to learn more about this unique system we obtained high resolution 
spectra during the eclipse of December 2001. We hoped to determine whether KH 15D is a 
radial velocity variable, whether there was evidence for additional light in the system beyond 
that of the K7 star, and whether the occultation had any effect on the spectrum. In fact, we 
found dramatic changes in the line profiles of Ha and H/5, as well as weak forbidden emission 
lines that become much more visible during eclipse. Evidently, the KH 15D system behaves 
like a "natural coronagraph" , allowing us to see details of its close circumstellar environment 
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during eclipse. Our evidence suggests that this WTTS is actively accreting gas and driving a 
bipolar outflow, although probably not at the rate of a typical CTTS. This calls into question 
the common practice of associating WTTS characteristics with the absence of an accretion 
disk. 



2. Observations and the Absorption Spectrum 

High resolution echelle spectra of KH 15D were obtained on the nights of UT 2001 
Nov. 29, when it was in its bright state just prior to eclipse, on UT 2001 Dec. 14, just 
past mid-eclipse, and again on UT 2001 Dec. 20 during egress (see Figure 1). These data 
were collected with the Ultraviolet- Visual Echelle Spectrograph (UVES) on the European 
Southern Observatory's Very Large Telescope (VLT), at Mount Paranal, Chile. The wave- 
length range is ~4780 to 6810 A. A 50 A gap centered on 5800 A is present due to use of 
the red arm which employs a mosaic of two 4096 x 2048 CCDs (D'Odorico et al. 2000). 
With a 1" slit, the spectral resolution is ~44,000. The spectra presented here were re- 
duced via the UVES pipeline. Information regarding the reduction process can be found at 
http://www.eso.org/instruments/uves/. As a check on this, the spectra were also reduced in 
the manner described by Valenti (1994). Both procedures make use of a sky subtraction al- 
gorithm. No significant differences were found between the data reduced by these techniques 
and the UVES pipeline reduction is adopted here. 

The UVES spectra confirm that, out of eclipse, the primary light source in the KH 15D 
system is a K7 WTTS. The 29 Nov. spectrum was visually compared in detail to 61 Cyg B, a 
typical K7 V, to look for any evidence of a second source of light in the system. No additional 
features were found. This spectrum was also compared to a rotationally broadened, synthetic 
K7 V template to determine its v sin(i), which we estimate to be < 5 km/s. The EW of the 
Lil 6707 feature is 0.401 A ± 0.001, consistent with the value of 0.47 A ± 0.05 measured on 
a low resolution spectrum by Hamilton et al. (2001). 

The out-of-eclipse spectrum was cross-correlated against HD 55999, a standard ob- 
served on the same nights, to determine the radial velocity of KH 15D. This was done with 
the IRAF 3 task FXCOR. Since the spectra were obtained on two CCDs covering different 
wavelength regions, a cross-correlation was performed on each, avoiding the emission line 
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regions of Ha, and H/3. A heliocentric radial velocity of +9.0 km/s ± 0.2, was found by this 
procedure for KH 15D on 29 Nov 2001. 

The egress spectrum was cross-correlated against the out-of-eclipse spectrum in the same 
manner to look for any radial velocity variation between 29 Nov. and 20 Dec. A difference 
in heliocentric radial velocity of 3.3 ± 0.6 km/s was found. The radial velocity of the star 
on 20 Dec. was, therefore, +12.3 km/s ± 0.6. Whether this detection of radial velocity 
variation means that the K7 star is a spectroscopic binary remains to be seen. Photospheric 
line profile variations are expected in a partially eclipsed star and the likely importance of 
scattered radiation near minimum light further complicates the interpretation. Clearly a 
more extensive radial velocity study of the system is needed and is underway (J. Johnson & 
G. Marcy, private communication). 

3. The Emission Line Spectrum 

In Figures 2-4, we show the Ha, H/3, and [01] 6300 emission-line profiles for KH 15D. 
Each spectrum has been flux calibrated to a relative scale using an R magnitude for the 
date, as given in Table 1. The quoted uncertainties on the magnitudes reflect a small degree 
of non-simultaneity in the spectral and photometric data as well as the need to transform 
from I (where the data were more numerous) to R. Since spectral and photometric data were 
taken within an hour of each other, except at maximum light (when the brightness is nearly 
constant with time) and since the color variation is quite small (~0.1 mag in R-I over the full 
brightness range), the uncertainty in R is only ~0.1 mag. Julian dates for each observation, 
measured I magnitudes, derived R and V magnitudes, measured EWs, and derived Ha and 
[OI] fluxes are listed in Table 1. The EW measurements refer to the entire profile, both 
red- and blue-shifted components combined. On Figs. 2 and 3, the arrows indicate where 
HI nebular emission could have affected the profile, whereas on Fig. 4 the arrow indicates 
where [OI] A6300 from the night sky may not have been removed completely. 

Figure 2 shows that dramatic changes occurred in the Ha emission line profile of KH 15D 
during eclipse and egress. Out of eclipse, the star appears to be a WTTS with an EW(Ha) 
~2 A, but with a double-peaked profile, having a central absorption and faint, but clearly 
detectable broad wings. A comparison with the Ha profiles of other WTTSs (Hartmann 
1982; Mundt et al. 1983; Finkenzeller & Basri 1987; Edwards et al. 1994; Reipurth et al. 
1996) shows that only 3 out of 19 such stars exhibited double-peaked profiles similar to KH 
15D, with UX Tau A being the most similar (see Reipurth et al. (1996)). Most WTTSs 
show narrow, single-peaked emission lines. It is also interesting that the central absorption 
feature appears to extend below the stellar continuum. This is unusual for any TTS, weak 
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or classical. During eclipse and egress, the "natural coronagraph" effect is clear. Near mid- 
eclipse, the EW of Ha grows to ~40 A, while the relative flux drops by ~50%. During egress 
the EW rises to ~70 A, as the flux increases, exceeding even its out-of-eclipse value. In 
all double-peaked profile is observed. The emission line profile during mid-eclipse 

extends to ± 200 km/s with significantly less flux in the extended wings, as compared to the 
profile during egress, which extends to ± 300 km/s. 

Figure 3 shows the corresponding H/3 line profiles. Again, dramatic EW and line profile 
changes are evident. Although the out-of-eclipse H(3 line profile is heavily affected by the 
underlying stellar absorption spectrum, it is similar to Ha. In the egress spectrum, however, 
the emission occurs primarily on the blue side with a wing extending to -300 km/s. The 
asymmetry in this line is much more striking than Ha. During mid-eclipse, although little 
H/? flux is present, the shape of the line profile appears similar to egress. The weak features 
visible in this profile at - 25 km/s are probably due to an improper background subtraction 
of the R/3 emission from NGC 2264. 

Figure 4 shows the emission-line profiles for the [OI] 6300 A line. The [OI] line, which 
is weak, but clearly discernable at maximum light, becomes prominent at mid-eclipse and 
during egress. The line flux seems to be about the same out of eclipse and during egress, 
although slightly higher at mid-eclipse. We caution against any extreme interpretation of 
this measurement. The profile during mid-eclipse was disturbed by an improper background 
subtraction, and we feel that our errors are most likely underestimated. However, this 
behavior would indicate that none of the [OI] emitting zone suffers variable occultation at 
the phases of our observations. This is fully expected given the spatial extent of the (bipolar) 
forbidden line emitting regions in CTTSs. The peak of the high velocity component of the 
[OI] emission in typical CTTSs originates at about 30 AU from the star (Hirth et al. 1997). 

The EW of the [OI] A6300 line in the out-of-eclipse spectrum is about 0.17 A, which is 
quite large for a WTTS. Only one of the ten WTTSs in Table 3 of Hartigan et al. (1995) 
has a detectable [OI] A6300 line (with an EW = 0.5 A). The others have upper limits of 
about 0.06 A. Most CTTSs in that table, on the other hand, have EWs of 0.5 - 3 A. The 
profiles obtained during mid-eclipse and egress suggest that we have two emission peaks, 
at about -20 and +18 km/s with emission wings extending -60 to +50 km/s, respectively. 
These profiles are quite different from the [OI] A6300 profile seen in most strong-emission 
CTTSs (with strong veiling), which often have a high- velocity component at -100 to -150 
km/s (resulting from the jet) and a low-velocity component at about -20 km/s (Hartigan 
et al. 1995; Hirth et al. 1997). In addition, the profile is quite different from that of a 
CTTS with small veiling and IR excess (see the bottom profile in Figure 11 of Hartigan et al. 
(1995)). These latter profiles are usually single-peaked and unshifted in velocity. Although 
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the profile shape for KH 15D is not fully clear due to imperfect subtraction of the emission 
of [01] at 6300 A in the night sky, it is most likely double-peaked. Such a profile can be 
most easily explained by a bipolar jet moving nearly perpendicular to the line of sight, quite 
consistent with our expectation that the disk associated with KH 15D is viewed nearly edge 
on, resulting in a very small radial velocity separation between the two jets. The profile 
catalog of Hartigan et al. (1995) contains several examples of bipolar jet sources (RW Aur, 
AS 353 A, DD Tau), where the two jet components are clearly separated due to favorable 
inclination angles. 

4. Discussion 

It is important to keep in mind that we do not yet know the geometry of the KH 
15D system. In particular, it is not certain whether the occultation proceeds along a line 
perpendicular to or parallel to the rotation axis of the K7 star (or neither) or to what extent 
the rotation axis, the presumed stellar magnetic axis and the orbital plane are aligned. 
It is not even known whether, relative to the system's center of mass, it is the K7 star 
or the occulting matter which is primarily in motion. Detailed modeling of the system 
must obviously await clarification of these basic issues. However, a preliminary qualitative 
interpretation of the emission line variations is possible and provided in this section. 

We believe that the observed spectral variations of KH 15D can be understood qualita- 
tively in terms of a weakly accreting, "scaled down" classical T Tauri star (see e.g., Muzerolle 
et al. (2001)) whose photosphere and magnetosphere are periodically occulted by a rela- 
tively sharp knife-edge (as described in Herbst et al. (2002), their Fig. 5). We expect such 
a star to have a bipolar jet region, which we assume, by analogy with the CTTS, is revealed 
by the forbidden emission line radiation. Assuming a jet velocity of 200 km/s (which is the 
average jet velocity for CTTSs derived by Hirth et al. (1994)) and adopting a radial velocity 
of ± 20 km/s for the [OI] A6300 peaks, we derive an inclination angle for the jet to the line 
of sight of 84°. This is consistent with a general picture for the system in which we view 
the K7 star close to the orbital plane of its circumstellar (or circumbinary) disk and the jets 
emerge roughly perpendicular to the disk plane, as for the star associated with HH 30 and 
other examples imaged by the Hubble Space Telescope (Ray et al. 1996). The absence of 
any significant variation in the profiles or flux of the forbidden line radiation during eclipse 
is consistent with the expectation, based on the CTTS analogy, that it arises at distances of 
tens of AU's from the star, beyond the region variably occulted. 

The behavior of the hydrogen lines is complex because some components do arise close to 
the star and, therefore, suffer variable occultation effects along with the photosphere. Since 
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we expect that the Ha line has a much higher optical depth than the B./3 line, it is obvious 
that any emission line region will be more extended in Ha than H/3. If we assume that the 
magnetic axis of the K7 star is tilted toward us at ~ 5-10°, which is a reasonable assumption 
supported by the [01] line profiles, at mid-eclipse, one can qualitatively understand the Ha 
line profile as resulting from low velocity material in the outer, more extended Ha emission 
region while the star and the 11(3 line forming region are obscured by the occulting disk 
material. This would also explain why there is almost no B./3 flux during mid-eclipse. During 
egress, we expect most of the Ha emission line region, as well as some of the H./3 emission 
line region, which is closer to the star, to be visible. The Ha line profile during egress has a 
low velocity emission peak with a central absorption feature similar to what is seen in the 
out-of-eclipse profile. Additionally, two "shoulders" appear along the profile at about ± 150 
km/s extending out to ± 300 km/s. These "shoulders" could be due to material rotating in 
the outer parts of the magnetosphere. However, given that the v sin(i) is measured to be < 5 
km/s, and that the magnetosphere only extends out to about 5-6 stellar radii, the maximum 
rotational velocity is about 5-6 v sin(i) or 25-30 km/s, making it difficult to explain the Ha 
emission-line profile with a rotating magnetopshere. 

A more attractive hypothesis is that the high velocity "shoulders" on the Ha line, 
so prominent in the egress spectrum, arise from material falling along magnetic accretion 
columns. This interpretation can also qualitatively account for the H/3 emission line profile 
during egress. Adopting values for the mass and radius of the K7 star from Table 1 of 
Hamilton et al. (2001), a free-fall velocity of about 380 km/s can be associated with 
material at the surface of the star. Since H/3 is produced much closer to the star, the blue- 
shifted wing extending to nearly -350 km/s could be representative of material accreting 
along magnetic field lines near the pole. The asymmetry seen in the B./3 emission line profile 
is most likely due to the fact that the star is slightly inclined toward our line of sight. The 
reappearance of both a blue and red wing to the Ha line during the early part of egress, when 
most of the stellar photosphere is still occulted, shows that both red-shifted and blue-shifted 
gas is present along the same line of sight towards the small portion of the photosphere and 
magnetosphere that is being uncovered first. This also supports an accretion interpretation 
as opposed, say, to a rotation interpretation for the high velocities. 



5. Conclusions 

It appears that KH 15D, although clearly a WTTS by the usual criteria of weak Ha 
emission, absence of UV or IR excess emission, and relative photometric stability, is still 
undergoing active accretion and driving a bipolar outflow. It provides a cautionary example 
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against assuming that all stars with WTTS characteristics no longer have accretion disks. 
The unique geometry of this system, in which a relatively sharp-edged occulting mask crosses 
the photosphere has created a "natural cornagraph" which enhances the visibility of the star's 
magnetosphere during eclipse. The occultation also evidently crosses the inner portion of 
the star's magnetosphere, where high velocity gas motions arise, probably from magnetically 
channeled accretion. Synoptic studies of this star may ultimately allow us to reconstruct 
aspects of the structure of its magnetosphere with spatial resolution that will be unobtainable 
in other objects for decades to come. 

We thank U. Bastian, M. Ibraghimov, J. Johnson, G. Marcy, F. Vrba, and J. Aufdenberg 
for helpful conversations and useful data on this star. W.H. acknowledges support from 
NASA through its Origins of Solar Systems program. 
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Table 1. Journal of Observations and Measurements. 



JD (2452...) I V(mag) R(mag) I(mag) EW {Ha) (A) EW (OJ) (A) Flux (Hc0 (x 10" 6 ) a Flux [OJ] (x 10" 6 ) a 


242.7296 
257.7840 
263.6871 


16.10 ± 0.02 
19.6 ± 0.1 
19.1 ± 0.1 


15.27 ± 0.02 
18.9 ± 0.1 
18.3 ± 0.1 


14.49 ± 0.02 
18.2 ± 0.1 
17.6 ± 0.1 


2.11 ± 0.08 
39.3 ± 1.5 
70.8 ± 5.1 


0.22 ± 0.04 
8.4 b ± 0.2 
3.1 ± 0.2 


1.6 ± 0.06 
1.1 ± 0.04 
3.4 ± 0.24 


0.17 ± 0.03 
0.24 ± 0.02 
0.15 ± 0.02 



a Calculated as EW (Ho) x 10-°- 4 * fl 

b This EW was measured with the spike removed. 
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Fig. 1. — The December 2001 eclipse of KH 15D; photometry from Herbst et al. (2002). 
Arrows indicate epochs at which the UVES spectra were obtained. 
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Fig. 2.— Ha profiles of KH 15D obtained with the VLT and UVES during the December 
2001 eclipse. The arrow indicates where the nebular emission of Ha is seen in the background. 
The velocities are shown in the reference frame of the star. 
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Fig. 3.— H/3 profiles of KH 15D obtained with the VLT and UVES during the December 
2001 eclipse. The arrow indicates where the nebular emission of H/3 is seen in the background. 
The velocities are shown in the reference frame of the star. 
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Fig. 4.— [01] 6300 A profiles of KH 15D obtained with the VLT and UVES during the 
December 2001 eclipse. The arrow indicates where the night sky emission of [OI] has not 
been removed completely. The velocities are shown in the reference frame of the star. 



